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ABSTRACT 

We estimate the enclosed mass profile in the central 10 pc of the Milky Way by analyzing the infrared 
photometry and the velocity observations of dynamically relaxed stellar population in the Galactic 
center. HST/NICMOS and Gemini Adaptive Optics images in the archive are used to obtain the 
number density profile, and proper motion and radial velocity data were compiled from the literature 
to find the velocity dispersion profile assuming a spherical symmetry and velocity isotropy. From these 
data, we calculate the enclosed mass and density profiles in the central 10 pc of the Galaxy using 
the Jeans equation. Our improved estimates can better describe the exact evolution of the molecular 
clouds and star clusters falling down to the Galactic center, and constrain the star formation history 
of the inner part of the Galaxy. 

Subject headings: Galaxy: center — Galaxy: structure — Galaxy: kinematics and dynamics 



1. INTRODUCTION 

The center of the Milky Way is the closest galactic nu- 
cleus, at a distance of ~ 8 kpc from the Sun (Ghez et 
al. 2008; Gillcssen et al. 2009), and thus is a good lab- 
oratory to study galactic nuclei with. Nonetheless, the 
Galactic center (GC) had not been studied well enough 
until 1990s, owing to large interstellar extinction between 
the GC and the Sun (Ay jej 30 mag; Rieke, Rieke, & Paul 
1989), and the limit in the near-infrared (IR) observing 
technology in the past. Advances in high- resolution near- 
IR instruments during the last two decades have yielded 
a wealth of information on the detailed structure of the 
GC. 

The most central part of the GC harbors a compact 
massive object, probably a super massive black hole 
(SMBH) with a mass of - 4 x 10 6 M (Ghez et al. 2008; 
Gillessen et al. 2009). Faint (compared to the stars dis- 
cussed right below) blue stars known as the "S-stars" 
or "S-cluster" are observed in the immediate vicinity 
(within 0.04 pc) of the SMBH. Krabbe et al. (1995) 
identified them as massive main-sequence stars with a 
spectral type of B0-B9. Further from the center, be- 
tween ~ 0.04 and ~ 0.4 pc, a few tens of OB supergiants, 
giants, and main-sequence stars are observed along with 
a pool of faint red stars. Paumard et al. (2006) argue 
that these young stars form two disk-like orbital config- 
urations, highly inclined and rotating counter directions 
to each other, but note also that Lu et al. (2009) find 
only one disk. 

Kinematical properties of old GC population have been 
studied as well, although not as intensively as for the 
young population. Several studies (e.g., Genzel et al. 
1996, 2000; Figer et al. 2003) presented and analyzed 
proper motion and radial velocity observations of the old 
stellar population. This old population of stars is well 
relaxed under the influence of SMBH and thus contains 
some information on the mass distribution and dynami- 
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cal environment around the SMBH. 

The presence of very young (< 10 yr) stars in the cen- 
tral parsec has been a puzzle since the strong tidal forces 
and magnetic fields in the GC as well as the elevated tem- 
peratures in molecular clouds form hostile star formation 
environment there. Inward migration of a star cluster 
that is formed far outside the central parsec through dy- 
namical friction has been proposed (Gerhard 2001) to 
solve this youth paradox, but it was shown that this 
scenario requires unrealistically extreme initial cluster 
conditions to explain the observed distribution of young 
stars in the central parsec (Kim & Morris 2003; Kim, 
Figer, & Morris 2004). Thus "in situ" star formation 
is a more likely scenario, and star formation through a 
gravitationally unstable gaseous disk around the SMBH 
(Nayakshin, Cuadra, & Springel 2007, among others) ap- 
pears to be the most promising model currently. 

The gas material that formed the young stars inside 
the central parsec had probably come from the farther 
galactocentric distances, and a ring of dense molecular 
gas extending 2-7 pc from the SMBH, called a "circum- 
nuclear disk" (CND; Christopher et al. 2005), is a good 
candidate for its origin. The CND itself shows some ev- 
idences of star formation (Yusef-Zadeh et al. 2008) as 
well. 

The initial configuration (e.g., radial and vertical sizes) 
of the star-forming gaseous disk in the central parsec will 
depend on the trajectory of the gas material infalling 
from the CND, which in turn depends on the shape of 
the gravitational potential in that region (i.e., the central 
few parsecs). The exact shape of the potential in the 
central few parsecs will also determine the degree of tidal 
compression and distortion during the infall of the gas 
toward the central parsec, which in turn will determine 
the star formation efficiency in the central parsec. 

The masses of the SMBH and its immediate vicinity 
(< 0.1 pc) have been estimated by analyzing the line-of- 
sight (LOS) velocities and proper motions of the S-stars 
(Ghez et al. 2008; Gillessen et al. 2009, among others), 



2 



OH, KIM & FIGER 



while the enclosed mass profile (EMP) between 10 and 
100 pc from the SMBH has been calculated by analyz- 
ing the LOS velocities of the OH/IR stars (Lindqvist, 
Habing, & Winnberg 1992). The EMP between - 0.1 
and ~ 5 pc has been studied either by interpreting the 
velocities of the CND as those of a rotating, circular ring 
(Guesten et al. 1987; Jackson et al. 1993; Christopher 
et al. 2005) or by using the LOS velocity dispersion of 
old, relaxed stellar population (Genzel et al. 1996, 2000; 
Figer et al. 2003; Schodel et al. 2007). The estimate 
from the former sensitively depends on the assumption 
of a rotating, circular ring, and all of the latter studies are 
based on the observations only out to ~ 0.8 pc (Genzel 
et al. 2000 make use of the number count measurement 
out to ~ 5 pc with the SHARP speckle camera on the 
3.5-m New Technology Telescope [NTT] obtained from 
the diploma thesis of Schmitt 1995, but the quality and 
reliability of this measurement is difficult to be assessed) . 
Thus the EMP estimate in the central few parsecs regime 
is still rather uncertain. 

In the present paper, we estimate the EMP in this im- 
portant region, between ~ 0.1 and ~ 10 pc from the GC, 
by analyzing the Hubble Space Telescope (HST) infrared 
photometry inside the central 5 pc along with the radial 
velocities and proper motions of the old stellar popula- 
tion in the same region. 

2. THE DATA 

The Jeans equation for a spherically symmetric sys- 
tem will be used in the present study to estimate the 
EMP in the GC, and for this, one needs number density 
and velocity dispersion profiles of a relaxed population. 
The former can be obtained from the stellar photometry, 
and the latter from the proper motion observations and 
infrared spectroscopy. 

2.1. Stellar Photometry 

Five near-IR images toward the GC taken with the 
NICMOS camera 2 (NIC2) onboard the HST, which are 
available from the HST archive, have been analyzed to 
obtain the stellar number density profile. We adopted 
F160W and F222M (similar to Johnson H and K) filter 
images observed in October 1997 and September 2002. 
Table Q] lists those image frames, and Figure [1] shows 
their sizes and locations, which cover the central ~ 4 pc 
of the GC (100"=3.88 pc at the assumed GC distance 
of 8 kpc). The five frames are roughly aligned on the 
Galactic plane, and the two larger frames are mosaiced 
ones each composed of 4 images. The pixel scale and the 
field of view of each NIC2 image (256 x 256 pixels) are 
0.076"/pixel and 19.2" x 19.2", relatively. 

We performed PSF photometry using the DAOPHOT 
package (Stetson 1987) within the Image Reduction and 
Analysis Facility (IRAF). 7 to 10 bright and relatively 
isolated stars were used to construct a point-spread func- 
tion (PSF) of each image. NICMOS PSFs have quite 
prominent secondary diffraction rings and radial spikes, 
and automated star-finding algorithms often identify the 
overlaps of two PSFs as stars. We removed these bogus 
stars by hand. We assume that the stars have the same 
intrinsic color (mFwow — itif222m)o — 0.25 mag (see 
Kim et al. 2005) and calculate the reddening of each 
star following an extinction law by Rieke et al. (1989). 



The average extinction in our F222M image is found to 
be ~ 3.2 mag. 

We carried out a completeness test by adding artificial 
stars to the observed images. We find that the 80 % 
completeness limit of the innermost F222M image is at 14 
mag, and use stars brighter than 14 mag for our analyses 
(the recovery fractions of the F222M images other than 
the innermost one is greater than 90 % at 14 mag). 

Estimation of an EMP requires a density profile of a 
tracer population that is dynamically relaxed. Most of 
the stars outside the central 0.4 pc of the GC are inter- 
mediate to old populations, while the stars inside 0.4 pc 
are a mixture of young and older populations (Krabbe 
et al. 1995). To cull out the young population from the 
older, we use the CO line strengths from the Adaptive 
Optics Demonstration Science Data Set of the Gemini 
telescope. 2 The Gemini GC demo data were observed in 
July and August 2000 with H, K' , A"-continuum (cen- 
tered at 2.26 fim with a bandpass of 60 A), and CO (cen- 
tered at 2.29 fxm with a bandpass of 20 A) filters. We 
performed PSF photometry for images 1 (roughly cen- 
tered at the Sgr A*) and 2 (centered about 20"away from 
the Sgr A*; each image has a field of view of 20" x 20"). 
We have calibrated the photometry of the two images 
using the stars that appear on both images, but they are 
not absolutely calibrated. 

Figure [2] shows K cont — CO vs. K cont diagrams for stars 
inside and outside the 0.4 pc radius from the Sgr A*. 
Outside 0.4 pc, the color-magnitude diagram forms a rel- 
atively narrow stream at the bright end (K cont < 16), 
but inside 0.4 pc, a separate population with redder 
K con t— CO colors (i.e., smaller CO absorption strengths) 
is seen. The small CO strength is an indication of young 
ages (< 10 7 yr), and we identify the stars in the box of 
Figure [5^, as the young population. We removed the 
young stars in the central 0.4 pc from our NICMOS 
photometry by cross-identifying the Gemini photometry 
against the NICMOS (the number of young stars was 34 
out of 1535 stars with F222M < 14 in our NICMOS pho- 
tometry). This way, we were able to have a collection of 
mostly intermediate to old stars that are brighter than 
14 mag in K' in the central 4 pc of the GC, and this 
sample will be used to produce the density profile of our 
tracer population. 

2.2. Stellar Velocities 

Genzel et al. (2000) compiled a homogenized data 
set of stellar velocities within the central 0.8 pc of the 
GC by combining various proper motion and LOS ve- 
locity data from the New Technology Telescope (NTT), 
the Keck telescope, and the MPG/ESO telescope on La 
Silla. Figer et al. (2003) reported LOS velocities of 85 
cool stars in the central parsec of the GC obtained with 
the Keck telescope. 3 We combined these two data sets to 
create a larger velocity sample of late-type (old) stars in 
the central parsec. 49 stars appear on both data sets, and 
we adopted the radial velocities from the Keck for those 



2 Available at http:/ /gemini.conicyt.cl/sciops/data/release_doc/manual.html 

3 Zhu et al. (2008) report the second epoch observations of the 
LOS velocities toward the same region with the same telescope, 
which mainly target the accelerations of the stars. The velocity 
data from the second epoch are very similar to the first epoch, i.e., 
Figer et al. (2003). 
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common stars. The total numbers of late-type stars in 
this sample are 80 for the proper motion data and 236 
for the LOS velocity data. 

Although not as much as found for young stars, old 
stars in the GC show some figure rotation as well. Fol- 
lowing Genzel et al. (1996), we have subtracted a rota- 
tional velocity of 



« roM = 24(AZ/5") u ' 4 km/s 



(1) 



from our LOS velocities (Al is the Galactic longitude 
offset from Sgr A*). 

The stars in this data set extend out only to ~ 0.8 pc 
from the Sgr A*, we add to this data set two LOS velocity 
dispersion values at ~ 1.3 pc and ~ 4 pc that Genzel et 
al. (1996) have obtained from the literature (see their 
Table 6). 

3. NUMBER DENSITY PROFILES 

Our photometry gives a distribution of projected stel- 
lar distances from Sgr A* (R) for the tracer population in 
the central 4 pc region, and the surface number density 
profile (E[i?]) obtained from this distribution is shown 
with asterisks in Figure [3l To obtain a spatial number 
density profile (n[r]; r is the spatial radius from Sgr A*) 
from S(_R), we model n(r) with three functional forms 
adopted in previous studies. 

Genzel et al. (1996) implemented the following func- 
tional form for the spatial density: 



n(r) 



So 



ro ! + (£)«• 



(2) 



Here, Sq, Vq, and a are the parameters to be found. 
The information directly available from observations is 
S(i?), and n(r) can be transformed into S(-R) by the 
Abel integral equation, 



E(i?) 



n(r) 



R 2 



: r dr. 



(3) 



We numerically integrate this equation when fitting our 
E(i?) data from the observations. Genzel et al. (1996) 
fixed a to be an asymptotic value at large r, 1.8, but 
here we keep it as a free parameter. This density model 
will be designated as G96. 

Saha, Bicknell, & McGregor (1996) modeled the sur- 
face density with 



S(i?) = So ( 1 + -2 



R 2 
R] 



(4) 



where So and Rq are to be fitted. This model is a variant 
of the modified Hubble law for elliptical galaxies, and has 
an analytical counterpart for the spatial density of 



n(r) 



R B(l/2,a) 



r 2 \ -(«*+l/2) 

1 + 



(5) 



where B is the beta function. Saha et al. (1996) con- 
strained a to be 0.4, but we leave it as a free parameter. 
This model will be called S96. 

Alexander (1999) adopted a broken power-law spatial 
density model: 



n(r) 



n (r/r ) a r < r 



no (r/r )' 



r > r , (3 > 1 ' 



(6) 



where no, ro, a, and (3 are the parameters to be fitted. 
Alexander (1999) assumed r = 0.4 pc (10") and /? = 1.8 
for his fits, but again, we keep them as free parameters. 
We call this model A99. 

When fitting the above functions to our NICMOS data, 
we try both x 2 test and Kolmogorov-Smirnov (KS) test 
(see, e.g., Press et al. 1992 for the latter), x 2 test is 
widely used as a standard statistical method for com- 
paring two distributions, but its result becomes sensitive 
on the choice of the number and ranges of the bins if the 
number of incidences in some of the bins is too small (say, 
less than ~ 10). Since the number of intermediate to old 
stars in our data sample is small at the very vicinity of 
the GC (i? < 5"), this may make the x 2 test somewhat 
unreliable. For this reason, we also use the KS test as a 
supplementary test. KS test utilizes the cumulative dis- 
tribution function instead of the histogram, and does not 
suffer the arbitrariness problem. Note that, however, KS 
test has its own shortcoming, and this will be discussed 
shortly. 

Figure [3] shows our best x 2 and KS fits of the above 
density models to NICMOS data (Table [21 lists our best- 
fit density model parameters). The % 2 fits show a good 
agreement between different density models although the 
innermost bin appears to be slightly over-fitted. The 
probabilities that the model has a different distribution 
from the observation are less than 2 % for all three mod- 
els, so the slight over-fit in the innermost bin is statisti- 
cally not important. 

The KS test results in similar fits for models S96 and 
A99, but it gives a rather large discrepancy at the inner- 
most bin for model G96. This is probably because (nor- 
malized) cumulative distribution functions always start 
with and end with 1, making the KS test rather in- 
sensitive at both ends of the distribution. Model G96 is 
the least flexible function near the core radius (ro), so 
it finds a bit difficult to adapt itself to a sudden change 
near the core radius seen in the observation. Nonethe- 
less, all three KS fits are consistent with the observation 
by better than 98 % just as in x 2 tests. 

We do not choose the best density model out of these 
six fits at this point. Instead, we will see how much 
difference is made to the final EMP estimates from these 
six models. 

Note that all of our best-fits give n oc r~ 15 or similar 
relations at large r. This is somewhat shallower than 
those by Schodel et al. (2007; r -175 ) and by Genzel et 
al. (2000; r -1 - 8 ). These previous estimates are based on 
the photometry with less sensitivity and/or smaller radii 
covered than in the present study, and this is thought to 
be the cause of such differences. 

4. VELOCITY DISPERSION PROFILE 

The radii of the stars with proper motion data in our 
sample range from ~ 0.05 to ~ 0.5 pc, and those with the 
LOS velocity data range from ~ 0.1 to ~ 5 pc. Thus our 
projected tangential and projected radial velocity disper- 
sions (<jt and (Jr, respectively) from the proper motion 
data only cover the deepest region, and our LOS velocity 
dispersion (a z ) covers much wider region except the in- 
nermost area. These dispersion profiles are shown with 
three different symbols in Figure 3J 

For a functional form of spatial (i.e., not projected) 
velocity dispersion (cr„), we adopt the following parame- 
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terization used by Genzel et al. (1996): 

<7v(r) 2 =<Tl+<T 2 (r/r 0v )-t. (7) 

For r Qv , we use the same ro or R of the density model 
that is used to fit the above equation to the observed 
velocity dispersions (see below). Thus only (Too, ao, and 
7 are the parameters to be found. 

We assume that the velocity distribution is isotropic. 
We do not try anisotropic velocity models in the present 
paper because the spatial coverages of our proper motion 
sample and LOS velocity sample overlap only marginally. 

Since we assume an isotropy for the velocity dispersion, 
the relation between the observed velocity dispersion and 
the spatial velocity dispersion is given by 

<M«)° - «dU« - = ^ f r m 

Thus obtaining a functional form for <r v (r) profile re- 
quires spatial and surface density profile information, 
and we try all of our six fits for E(i?) (and its corre- 
sponding n[r\) obtained in the previous section. Figure 
0] plots our model (eq. [5]) fits to the observed velocity dis- 
persions with three different density profiles (G96, S96, 
and A99) and two different statistical tests (% 2 and KS). 
All six fits result in very similar velocity dispersion pro- 
files (best-fit velocity dispersion parameters are listed in 
Table©. 

5. ENCLOSED MASS PROFILE IN THE GALACTIC 
CENTER 

The Jeans equation is the first moment of the collision- 
less Boltzmann equation, and it gives a relation between 
the enclosed mass of a system and the velocity disper- 
sions. For a spherically symmetric, rotating system with 
an isotropic velocity distribution, the equation becomes 

GM(r) _ a 2 f rfln[n(r)] d\n[a v {rf] \ 

— V ro t\r) -\-o v \r) < > , 

r L a in. 7" ainrj 

where M is the enclosed mass and v ro t is the rotational 
velocity. The latter as a function of r can be obtained 
from the observed, mean LOS velocities as a function of 
the longitude offset from Sgr A* by the following relation: 

v rottZ (Al): 

M _ r f°° d ( Z(Al)v rot , z (Al) \ dAl 

Vrot[r) nn(r) J r dAl { Al J ^AP~r^ 

(10) 

For the density profile in this Abel transform, we use the 
density model for old population in Genzel et al. (1996) 
instead of those obtained in the present study because it 
better represents the whole inner bulge (several tens to 
hundreds of parsecs) in which the figure rotation takes 
place. 

Figure [5^ shows the EMPs obtained from the above 
Jeans equation with our density and velocity dispersion 
fits for the intermediate to old stellar populations. These 
profiles are quite similar to each other, but quite differ- 
ent from the two previous studies, Genzel et al. (2000) 
and Schodel et al. (2007). Our EMP is similar to the 
former at larger radii but is considerably larger than the 
former. The latter has much smaller EMP in the larger 



radii because it is based on a rather simple approxima- 
tion that the LOS velocity dispersion in the outer region 
is constant. 

We do not plot our EMPs inside 0.2 pc. Our EMP 
estimates at this inner region are not reliable because 
there are not many bright, intermediate to old stars in 
this region (this "hole" of old stars is probably due to 
frequent close encounters with other stars; see Genzel et 
al. 1996). 

For a smooth convergence of the EMP to the mass of 
the SMBH, we extrapolate the EMP inward from the 
radius at M = 5 x 10 6 M Q assuming that the extended 
mass (the mass excluding the SMBH) distribution follows 
a power-law function and that the mass of the SMBH 
is 4 x 1O 6 M (Ghez et al. 2008; Gillessen et al. 2009). 
FigureOJ) (and TableS]) shows such inward extrapolation 
of our x 

best-fit of A99 density model (A99/x 2 )- We 
choose this model as our canonical fit as it describes the 
observed density profile the best and its EMP roughly 
represnts the average of all our 6 models. 

Local densities of the extended mass (the mass other 
than the SMBH) can be obtained by differentiating the 
EMP. The density profile for our A99/x 2 model is given 
in Figure [5] and Tabled] Again, our result is quite close 
to that of Genzel et al. at larger radii, but closer to that 
of Schodel et al. at smalle radii. 

6. SUMMARY AND DISCUSSION 

We have estimated the EMP in the central 10 pc of 
the Milky Way by analyzing the infrared photometry and 
the velocity observations of dynamically relaxed stellar 
population in the Galactic center. HST/NICMOS and 
Gemini images were used to obtain the number density 
profile of the relaxed population, and the LOS velocities 
and proper motion data of the same population were used 
to calculate the EMP from the Jeans equation assuming 
a spherical symmetry and velocity isotropy. The newly 
obtained EMP is larger than the previous studies at 0.1 < 
r < 10 pc, which is consistent with the most recent value 
for the mass of the SMBH. 

As discussed in §1, the enclosed mass and density pro- 
files in the central few parsecs will determine the exact 
morphological evolution of the molecular clouds and star 
clusters that are moving down to the GC. Our larger 
EMP implies that star clusters and molecular clouds in 
the central few parsecs will have smaller tidal radii and 
shorter orbital periods than previously expected. 

Degenerate stellar remnants such as stellar mass black 
holes and neutron stars are thought to be segregated 
in the GC due to their gradual dynamical friction with 
lighter field stars (Morris 1993). Our density profile can 
be used to constrain the amount of segregated degenerate 
objects. Since the current reservoir of degenerate stars 
is a result of continuous star formation in the GC (Figer 
et al. 2004), a detailed study on the stellar populations 
in the GC based on our EMP estimate will constrain the 
star formation history in the GC as well. 

The anonymous referee suggested density profile func- 
tions that have an inner radial cutoff to better describe 
the "hole" of old stars. We do not try such functions 
here because they will not significantly change our re- 
sults on the EMP and the nature of the hole is beyond 
the scope of the present paper, but we do believe that 
such functions can give some constraints on the size of 
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TABLE 1 

HST/NICMOS Datasets Analyzed in the Present Study 



Dataset Name 


Exposure (Sec) 


Observation Date 


F160W 






N49Z010H0 


25.9 


October 1997 


N49Z01030 


55.9 


October 1997 


N49Z01070 


55.9 


October 1997 


N49Z010B0 


55.9 


October 1997 


N49Z010F0 


55.9 


October 1997 


N49Z02030 


55.9 


October 1997 


N49Z02070 


55.9 


October 1997 


N49Z020B0 


55.9 


October 1997 


N49Z020F0 


55.9 


October 1997 


N6LO02050 


24.0 


September 2002 


N6LO03050 


24.0 


September 2002 


F222M 






N49Z010I0 


71.9 


October 1997 


N49Z01040 


207.9 


October 1997 


N49Z01080 


207.9 


October 1997 


N49Z010C0 


207.9 


October 1997 


N49Z010G0 


207.9 


October 1997 


N49Z02040 


207.9 


October 1997 


N49Z02080 


207.9 


October 1997 


N49Z020C0 


207.9 


October 1997 


N49Z020G0 


207.9 


October 1997 


N6LO02040 


71.9 


September 2002 


N6LO03040 


71.9 


September 2002 



TABLE 2 
Surface Number Density Fits 



G96 



S96 



A99 



G96 



S96 



A99 



Eu- 
ro 



0.11 
1.50 



Ro 

a 



Fit 
1130 
0.13 
0.25 



no 
ro 
a 





994.7 
0.17 
0.45 
1.47 



ro 



KS Fit 



0.44 
1.48 



Ro 
a 



1100 

0.14 
0.25 



n 
r 
a 




1168.1 
0.15 
0.59 
1.48 



Note. — So are in units of pc , no in pc , and ro & -Ro in pc 



TABLE 3 
Velocity Dispersion Profile Fits 



G96 



S96 



A99 



G96 



S96 



A99 



Coo 
°0 

7 



"3W" 

253 
0.77 



X 1 Fit 
31.3 
237 
0.77 



-287T 

209 
0.73 



170 
1.11 



KS Fit 
44.8 
201 
0.99 



210 
0.81 



Note. - 



(Too and ijo values are in units of km/s. 
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TABLE 4 

Enclosed Mass and Density Profiles in the Galactic Center 



r 


M 


P 


(PC) 


(M Q ) 


(M /pc 3 ) 


0.1 


4.15c+06 


2.49c+07 


0.2 


4.40c+06 


5.87e+06 


0.3 


4.73e+06 


3.25c+06 


0.5 


5.54e+06 


1.33e+06 


0.7 


6.33e+06 


6.30e+05 


1.0 


7.47e+06 


3.04e+05 


2.0 


1.16c+07 


8.88e+04 


3.0 


1.65e+07 


4.76e+04 


5.0 


2.90c+07 


2.29e+04 


7.0 


4.48c+07 


1.42c+04 
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100 




-100 i i i i i i i i i i i i i i i i i i i 

100 50 -50 -100 

RA offset from Sgr A* (") 

Fig. 1. — Locations and sizes of the HST/NICMOS images used in the present study. The cross at the center indicates the Sgr A*. The 
five frames are roughly aligned along the Galactic plane, and the two larger frames are mosaiced ones each composed of 4 images. 



12 



a) r<10" 



16 ■ 



18 ■ 



20 - 



22 

-2.0 -1.5 



Young 
• Population 



-1.0 -0.5 

K,„„,-C0 



0.0 



12 



0.5 



b) r> 1 0' 



14 - 



16 - 



20 - 



22 

-2.0 -1.5 



.0 -0.5 
K r „ nt -C0 



0.0 0.5 



Fig. 2. — K con t— CO vs. K con t diagrams for the stars in r < 10" (a) and r > 10" (b) regions from the Gemini Adaptive Optics data. The 
magnitudes presented here are not calibrated ones (we find that our Gemini K con t magnitudes arc ~ 3.6 mag fainter than the NICMOS 
F222M magnitudes on average). A separate population of stars is seen in the area defined by K con t — CO > —0.75 and K con t < 17.5 
(denoted with the straight lines) at r < 10". This population is mostly composed of young stars (see the text). 
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Fig. 3. — Best-fit surface density profiles for density models of Genzel et al. (1996; G96), Saha et al (1996; S96), and Alexander (1997; 
A99) using the \ 2 test (a) and the Kolmogorov-Smirnov test (b). The surface densities from our analysis of HST/NICMOS images (after 
the subtraction of young population) arc shown with asterisks and 1-a error bars. 




Fig. 4. — Best-fit velocity dispersion profiles for density models of Genzel et al. (1996; G96), Saha et al (1996; S96), and Alexander 
(1997; A99) using the \ 2 test (a) and the Kolmogorov-Smirnov test (b). The velocity dispersion data, which we obtained from Genzel et 
al. (1996, 2000) and Figer et al. (2003), are shown with symbols and 1-cr error bars. 
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0.1 1.0 10.0 0.1 1.0 10.0 

r (pc) r (pc) 

Fig. 5. — (a) Enclosed mass profiles from the Jeans equation for density models of Genzel et al. (1996; G96), Saha et al (1996; S96), 
and Alexander (1999; A99) and statistical tests of x 2 and Kolmogorov-Smirnov (KS) methods, (b) Enclosed mass profile of A99/x 2 model 
modified to converge to the SMBH mass of 4 X 10 6 Mq at r = pc (our best-fit). Also plotted are the enclosed mass profiles estimated by 
Schodel et al. (2007) and Genzel et al. (2000). 




0.1 1.0 10.0 

r (pc) 

Fig. 6. — Density profile of the modified A99/x 2 model (our best-fit; see Fig. [5p) along with those estimated by Schodel et al. (2007) 
and Genzel et al. (2000). 



